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We investigate the shear viscosity of thermalized photons in the quark gluon plasma (QGP) at 
weak coupling and A/” = 4 super Yang-Mills plasma (SYMP) at both strong and weak couplings. 

We find that the shear viscosity due to the photon-parton scattering up to the leading order of 
electromagnetic coupling is suppressed when the coupling of the QGP/SYMP is increased, which 
stems from the blue-shift of the thermal-photon spectrum at strong coupling. In addition, the shear 
viscosity rapidly increases near the deconfinement transition in a phenomenological model analogous 
to the QGP. 


The electromagnetic (EM) signatures such as leptons 
and photons produced in thermal plasmas play an im¬ 
portant role in probing the properties of the plasmas, 
which have been widely investigated in relativistic heavy 
ion collisions and cosmology. Peculiarly, recent observa¬ 
tions of large elliptic flow of direct photons comparable 
to that of hadrons in heavy ion collisions m lead to the 
tension between experimental results and theoretic pre¬ 
dictions Em]. This inconsistency also stimulates many 
studies for the electromagnetic probes in the QGP with 
distinct mechanisms [SHU]. Although the photons gen¬ 
erated in heavy ion collisions hardly reach thermal equi¬ 
librium due to short lifetime of the QGP, some of the 
photons or leptons produced in the cosmic plasma with 
longer lifetime may reach equilibrium. In such a case, e.g. 
in the quark epoch, the cosmic plasma could be regarded 
as a QED-bQGD plasma. 

It is found that the quantum electrodynamic (QED) 
plasma is rather viscous and the thermal leptons or pho¬ 
tons barely interact with each other after being emit¬ 
ted from the plasma dSHm, whereas the scenario in 
the QGP is not so obvious since the interaction be¬ 
tween the leptons/photons and the medium could be 
enhanced by the lepton-quark and photon-parton scat¬ 
terings. Particularly, the photon-parton scattering dom¬ 
inates the photon-photon/lepton interaction in the or¬ 
der of EM coupling. Moreover, in the strongly coupled 
QGP (sQGP), the non-perturbative effect in the color 
sector may modify the interaction with the EM sector. 
We thus consider the scenario when leptons/photons and 
the QGP reach thermal equilibrium with same temper¬ 
ature, while they may not share same transport coef¬ 
ficients due to different couplings in the EM and color 
sectors. Theoretically, it has been shown recently that 
the shear viscosity of thermalized leptons in the QGP 
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could be suppressed by the lepton-quark scattering com¬ 
pared with the case in the QED plasma [TH|. There¬ 
fore, it is tentative and imperative to analyze the photon 
transport in the QGP. In this paper, we follow the ap¬ 
proach in [TH] to compute the shear viscosity of photons 
stemming from the photon-parton scattering in the QGP. 
We apply the relativistic Boltzmann equation to describe 
the dynamics of thermalized photon, while the collisional 
terms could be obtained from both perturbative or non- 
perturbative approaches depending on the coupling of 
the color sector. However, due to the limitation of lat¬ 
tice simulations for real-time observables, we may resort 
to the gauge/gravity duality [TM^ . the correspondence 
between d-dimensional strongly coupled gauge theories 
in large A), and d 1-dimensional supergravity, to quali¬ 
tatively delineate the non-perturbative properties of the 
sQGP. One will see that the corresponding spectral func¬ 
tion beyond the hydrodynamic regime at strong coupling 
yields nontrivial suppression of the shear viscosity of pho¬ 
tons. 

Gonsidering only particle scattering, the relativistic 
Boltzmann equation of photons can be written as 

x) = -f{p, x)r> (p) + (1 + f{p, a;))r< (p), (1) 

pU 

where 

r<(»(p) = (2^)3^^. (2) 

Here represent the production (absorption) rates 

per spacetime volume and f{p,x) denotes the distri¬ 
bution function of photons in phase space. Based on 
a naive power counting, the electron-quark/electron- 
scattering cross section is O(e^) [52], whereas the photon- 
quark/gluon-scattering cross section is 0{e^). Intu¬ 
itively, one could expect more substantial influence from 
the sQGP on photon transport. To study the photon 
transport, we apply the linear response theory to in¬ 
troduce the local fluctuations of photons slightly away 
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from equilibrium in the vicinity of the rest frame but 
keep the sQGP in equilibrium. Our approach here is 
distinct from studies of the transport properties of a 
single system with multiple components such as weakly 
coupled QGP na [H], where the different quarks and 
gluons share the same form of fluctuations. In gen¬ 
eral, the fluctuations of photons also induce the back- 
reaction on the sQGP, while the back-reaction should be 
higher-order corrections in e. Moreover, by implement¬ 
ing the non-equilibrium fluctuation-dissipation theorem 
in holography [52], one finds that the fluctuations of the 
production and absorption rates in Q cancel each other 
[K5] . Based on the reasons above, we discard the fluctu¬ 
ations of the color sector. In other words, we treat pho¬ 
tons as probes; the photon transport here only encodes 
the equilibrium properties of the sQGP. By introducing 
/(p, x) = nb(p, x) + Sf{p, x) with nb{p, x) being the ther¬ 
mal distribution of bosons, we find 


-^df,f{p,x) = - -- i-Sf{p,x), 

pD ^ pV 


( 3 ) 


Alternatively, one may rewrite it in terms of the produc¬ 
tion rate. 


1 

“ 120T7r4 


d\p\\pfnt{p,xf{l + nb(p,x)) 



.( 8 ) 


In fact, following |23j . one can further derive a more gen¬ 
eral expression of the shear corrections, which takes the 
form 

Sf^Hk) = [ = 2G^{k)d^^u^\ (9) 

J (27r)-^p‘J 


where the bar corresponds to the Fourier transform with 
respect to x and k denotes the dual momentum. Here 
is related to the retarded correlation function of the 
energy stress tensor, which reads 
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d^P IpP _ npjl + Up) 

(27r)3 T AfcO_fc.^+ Mnj(P°)) 


410) 


where n^ (p°) denotes the trace of the light-like retarded 
EM-current correlator of the sQGP. To obtain the second 
equality above, we utilize the following relations. 


nb{p,x) dT^ 
d^p 1 -I- nb{p, x) d^p 


nb{p,x) 

(27r)3p0 


Im [n«(p°)] , (4) 


where the photon-production (absorption) rates are as¬ 
sociated with the retarded EM-current correlators. The 
Boltzmann equation ([^ implies 


'im [n«(p°)] 


= 2p\ip^) 


,0\-l 


( 5 ) 


where r-y denotes the relaxation time of photons, which is 
shown to be inversely proportional to the light-like spec¬ 
tral function x(p°). Note that the relaxation time here 
depends on the energy of photons. For the low-energy 
photons, one simply finds » l/(2crc), where CTc corre¬ 
sponds to the DG conductivity of the sQGP. 

Now, by taking the shear fluctuation as 
Sf{p,x) = {1 + Ubip, x))nb{p,x)^^^f^ (6) 


where dt^iUj) = {diUj + djUi) /2 and T denotes the equi¬ 
librium temperature and represents the collective ve¬ 
locity of photons. After solving ([^, we obtain 


B{p) = 


Im [n^(pO)] ’ 


( 7 ) 


which gives rise to the shear viscosity of photons. 


_ 1 /■,, ,\P?nb{p,x){l+nb{p,x)) 

60T7r2 7 'P' Im[n«(pO)] 


where v = p/\p\- The shear viscosity is defined as 
r]^ = G^(0), which agrees with the expression in 
According to [531 HD, one can further study the shear 
relaxation time = —i/{k°)i with (/c°)i being the first 
pole of G^{k^,k = 0). Although the shear relaxation 
time is important for understanding shear corrections, 
numerically finding the poles from (101 is more involved. 
We thus only focus on the shear viscosity in this paper. 


The general expression in Q is applicable for both the 
strongly/weakly coupled QGP or different types of media 
coupled to photons. We focus on the universal feature of 
rjj affected by the coupling of the medium. Thus, we will 
firstly investigate the in the Af = 4 super Yang-Mills 
(SYM) theory despite its difference from QGD, where 
the EM-current correlators of the SYM plasma in both 
weakly coupled and strongly couped scenarios have been 
studied [25]. In addition, the photon emission at finite 
but large t’Hooft coupling A in the SYM plasma was an¬ 
alyzed in holography EllIT], which allows us to further 
explore the change of r]^ in large A. The EM-current cor¬ 
relator in the weakly coupled QGD at finite temperature 
could be found in [HHI]. Nonetheless, to qualitatively 
capture the non-conformal effect near the critical tem¬ 
perature in the strongly coupled scenario, we will em¬ 
ploy a phenomenological model in holography to mimic 
the sQGP Hi HD. 


We firstly consider the photon transport in the Af = 4 
SYM plasma at strong coupling and large Nc limit, where 
the thermal-photon production rate was computed in 
[55] , where two of the Weyl fermions have electric charge 
±1/2 and two complex scalars have electric charge 1/2 
[54] . Although the charge assignment is not unique, the 
corresponding electric conductivity roughly matches the 
lattice simulation around T ~ S.ST^ [30] (see the compar¬ 
ison in m)- From [55] , the trace of the light-like spectral 
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function reads 


x{p°) 


- 2 a;"^Ini n^(p°) 


( 11 ) 


1 


771 / 2 — (l-|-i)cj 2-|-(l — -| • -| \ 

2F1 ( ^2 S i 1 “ “IJ 


where u) = / (27rT) and the hat corresponds to the nor¬ 
malization by e^iV^T^/S. By using ® and (11 1 and per¬ 


form the numerical integration of p, one finds 


„SYM 

h 


(A = 00 ) = 


i.46r3 

e'^m ■ 


( 12 ) 


In comparison with the shear viscosity of electrons rje 
in the strongly coupled SYM plasmas [E], the f]^ is 
C>(e^ ln(l/e)) suppressed [^. Since we now treat ther¬ 
mal photons as an independent system separated from 
the SYM plasma, it could be intriguing to compare the 
^SYM where denotes the entropy density of pho¬ 
tons, with the ratio to shear viscosity and entropy density 
of the strongly coupled SYM plasma itself, which is equal 
to l/( 47 r) and known as the Kovtun-Son- Starinets (KSS) 
bound [321 [33] • By simply using the entropy density of 
an ideal photon gas, = 47 r^T^/ 45 , we find 


= 00 ) 


25 
47r ’ 


(13) 


where we take Nc = 3 and Although the ratio 

is much larger than the one for the SYM plasma itself, 
it is highly suppressed compared with that in the QED 
plasma or in the weakly interacting case as we will see. 


Subsequently, we will evaluate {X) with large but 

finite A, which allows us to track the coupling depen¬ 
dence of . Considering the inclusion of full 0{a'^) 

type IIB string theory corrections, the strongly coupled 
SYM theory receives the 0{X~^^^) correction in hologra¬ 
phy. The thermal-photon production in the SYM plasma 
with such a correction was studied in |26j . Technically, 
the computation of the EM-current correlator requires 
solving a Schrodinger equation. 


dl'i’iu) 

V{u) 


V{u)'i’{u) 
1 


(14) 


fiu) 




- 11700 


-b2098482M2 _ 4752055 u^- b 1838319u® 
-ba;2u(-16470 -b 245442^^ -b 1011173u^) 


where u corresponds to the bulk direction with u = 0 
being the boundary and f{u) = 1 — \ with Uh being 

the position of the horizon in the dual geometry. Here 
7 = |C(3)A“^/^ denotes an expansion parameter with C, 
being the Riemann Zeta function. Eq. (141 has to be 


solved perturbatively in 7 with the in-falling boundary 
condition at the horizon; we thus write 

4'(u) = (1 - u)5(i-*‘")($o(w) + 7$i(u)). (15) 


Up to 0(1), solving (14) yields 


$o(w) = (l+u)5(i-“) 

"2 — (1 -b i)ui 2 -b (1 — i)oj 

I" 


X 27^1 


; 1 - jw; 


(16) 
1 — u' 


which is simply the solution giving rise to ( 11 ) at infinite 
coupling. The next-leading-order ( 0 ( 7 )) equation has 
to be solved numerically. Now the trace of the spectral 
function is given by 


X(P°) = 4Im 




4'n 


^'1 


4';, /265 




4', 


4'n 




4'r 


im 


u—O 


where the prime denotes the derivative with respect to 
u. As shown in |26j . the reduction of A results in the 
enhancement of the peak of the spectrum and the shift 
toward the infrared region. Although the DC conduc¬ 
tivity and the maximum amplitude of the spectrum are 
enhanced for smaller A jjSYM ijiQj. 0 ases when A de¬ 
creases as shown in Figj^ where the divergence of at 
A « 50 is due to the breakdown of the perturbative ex¬ 
pansion. On the contrary, by using the result in [18j . 
where /T^ is inversely proportional to the ratio of 

DC conductivity of the SYM plasma to temperature, the 
shear viscosity of electrons is suppressed when A 

decreases as shown in the same plot [SB]. In Figj^ we 
also present the shear viscosity at weak A, where the spec¬ 
tral functions and DC conductivity are obtained from the 
perturbative calculations in the SYM plasma [25] and in 
the QGP |31H5B] . Here we cut the results at A « 2.5, 
where the Uc of QGP starts to decrease when increas¬ 
ing A, which may imply the breakdown the leading-order 
computation [57]. Also, the Uc of the SYM plasma be¬ 
comes singular at A = 4.5. It turns out that the in¬ 
crease of with respect to the decrease of A is as 

well found in the weakly coupled scenario. Analogously, 
the decrease of with respect to the decrease of A 

is observed. Unlike in the SYM plasma, where the 

low-energy scattering dominates and thus the increase of 
the DC conductivity of the plasma reduces , the 

r^SYM jg governed by the full spectral function. It is in 
fact the blue-shift of the spectrum led by the increase 
of A found in both weakly coupled and strongly coupled 
scenarios causes the suppression of ij-y. Based on the qual¬ 
itative features of the shear viscosity at weak and strong 
couplings, we may expect monotonic decrease of rij and 
monotonic increase of r]e with respect to A. Note that 
^SVM (2 5 ) ~ { 00 ), which may further suggest 

that could almost saturate the value at A = 00 at 

intermediate A. 

In the following, we investigate the shear viscosity of 
photons in the sQGP via holographic QCD. Although 
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FIG. 1: The shear viscosity ratios of photons (red-solid) and 
leptons(blue-dashed) from holography in the strongly inter¬ 
acting SYM plasma with large t’Hooft coupling. 







FIG. 2: The shear viscosity ratios of photons (red-thick) 
and leptons (blue-dashed) for the weakly interacting SYM 
plasma. The shear viscosity of photons (orange-thin) and 
leptons (black-dotted) in the weakly coupled QGP are also 
scaled by those in the SYM plasma at A = oo. 


the N = A SYM theory and QCD may have similar fea¬ 
tures at intermediate temperature, the difference of two 
theories will become substantial near the critical temper¬ 
ature, where the conformal anomaly of QCD results in 
the deconfinement transition. In order to break the con¬ 
formal symmetry in holography, one has to incorporate 
the running of coupling, which could be characterized by 
a bulk scalar field with a corresponding potential in the 
gravity dual in the bottom-up approaches. This type of 
models may be regarded as an effective theory of QCD in 
the IR regime. Here we will employ the phenomenological 
model in Hi 133, which generates thermodynamic prop¬ 
erties similar to the results in lattice QCD. In [28l I37j . 
the gravitational action in the Einstein frame takes the 
form 






V{^) , 


(18) 


where G represents the determinant of the spacetime 
metric and V{(j)) denotes the scalar-field potential. 
Here = IGtt^L/TV^ with L being the AdS radius. The 
model is usually solved in the Gubser gauge with the 


following ansatz of the metric, 


ds^ = (-h{cl>)df + dx^) + (19) 

h{(p) 

where the scalar field </> is set as the fifth coordinate in 
the bulk. The blackening function h{(j)) should have a 
simple zero at the horizon (j) = (j)^ and the metric should 
recover AdS^ geometry near the boundary at (?!)—>■ 0. The 
temperature and entropy density are given by 

y ^ ^ g ^ ^ ( 20 ) 

Att Kj 

where the prime above denotes the derivative with re¬ 
spect to (j). In light of [29], we choose the following po¬ 
tential, 

V{(j)) = -12 cosh 7 ()) -I- 62 ')'^ + + becj)^, (21) 

where 7 = 0.606, &2 = 0.703, 64 = —0.12, feg = 0.0044, 
and the AdS radius L = 1. We will not present the details 
for solving the metric, which could be found in H5j. 


Now, we have to coupled an U(l) gauge field to the 
gravitational action. The relevant term which follows 
the linear response of the U(l) current can be written as 


Sf = 
where Fmn 


"2^ / 

denotes the field strength. We will take 
- sech(ai<))) 

/(*r) 2 ’ 


( 22 ) 


(23) 


which is introduced in HS] to fit the electric susceptibil¬ 
ity from lattice simulations for T < 1.5Tc |S3] , where 
ai = 0.4 and 557 is a dimensionless constant associated 
with the overall amplitude of the observables. More im¬ 
portantly, the choice qualitatively captures the increase 
of CTc near Tc and the saturation of the SYM result around 
T « 3.5Tc in comparison with the lattice simulations m- 
We then solve the Maxwell equation numerically with the 
incoming-wave condition at the horizon. From the stan¬ 
dard AdS/CFT prescription |35|, the trace of the spectral 
function reads 


p2 /V2 

xb") = 


lim f/(^)VAGG‘^'^G“^^ 
4>^o y Ei 


,(24) 


where Ei = |p|Ai. In Figj^ we illustrate the thermal- 
photon spectra at different temperature in strong cou¬ 
pling, where 




dT</dp° 

aEMN^T^' 


(25) 


Here we extract from the minimum of the speed of 
sound 133 and fix the overall amplitude of the spectrum 





























5 


(dr,/dp")N„„ 




FIG. 3: Solid curves from bottom to top correspond to the 
photoemission rates in holographic QCD in strong coupling 
at T =(1, 1.08, 1.26, 1.54, 3.57, 7.4) Tc, respectively. Dashed 
curves from bottom to top at p°/T = 2 correspond to the 
rates at weak coupling with as = 0.05 and as = 0.1. 


FIG. 4: The solid red curve corresponds to normalized con¬ 
ductivity obtained from the holographic QCD and the dashed 
line corresponds to the one from strongly coupled A/” = 4 SYM 
plasma. The blue points represent the lattice simulation for 
the light-fstrange case in |39| . Here Cem = ji. 


with a proper choice of by matching the result at 
T = S.SdTc, where the spectra at higher T start to sat¬ 
urate, and the one from the A/” = 4 SYM plasma. Such 
a choice yields the conductivity shown in Fig|^ in com¬ 
parison with the one in the SYM plasma [5S] at A = oo 
and the lattice simulation in [30112S]. InFigd we also 
show the results from the weakly coupled QCD at high 
temperature [MUSS], where we take Nf. = 3 and Nf = 3. 
Similar to the scenario in the A/” = 4 SYM plasma, one 
finds substantial blue-shift of the photon spectrum at 
strong coupling. In fact, the blue-shift may be fore¬ 
seen in the long-wavelength limit. Up to the leading 
order of the small-momentum expansion of the light¬ 
like spectral functions, the AdS/CFT calculations yields 
(x/T^) ^ p°/T, which follows the typical pattern of hy¬ 
drodynamics, whereas the perturbative QCD and Af = 4 
SYM theory lead to (x/T^) ^ ln(p°/T) and 

(x/T^) ~ ln(p°/r), respectively. 

However, to locate the peak of the spectrum beyond 
the hydrodynamic region, one has to perform the com¬ 
putations of full spectra. In the holographic model, we 
find that the peaks of the spectra at different temper¬ 
ature almost locate at the same momentum scaled by 
the temperature. Furthermore, the normalized spectra 
with high temperature far from Tc almost coincide. This 
finding agrees with the saturation of the ratio of con¬ 
ductivity to temperature found in the same model |2^ . 
After obtaining the photon spectrum, we can employ (1^ 
to evaluate the shear viscosity of thermal photons in the 
holographic QCD Since the shape of the spec¬ 
trum is nearly unchanged by varying T, scaled 

by rapidly increases near Tc as shown in Figj^ which 
is qualitatively analogous to the behavior of the shear 
viscosity of thermal leptons [T5| . 

In summary, we found two salient features of the pho¬ 
ton transport in the sQGP. First, albeit it is only shown 
in the A/" = 4 SYM plasma, thermal photons in the sQGP 
become more fluid-like when the coupling increases. Sec- 


Ir/e 


hQCD 



FIG. 5: The solid and dashed curves correspond to the ratio 
of 77 t and the ratio of rjc. 


ond, the shear viscosity of photons in the sQGP increases 
near the deconfinement transition. Nevertheless, both 
the shear viscosity and relaxation time of thermalized 
photons depend on the electric conductivity of the QGP. 
According to [30] , CTc in the QGP from recent lattice sim¬ 
ulations is not very far from the value of the strongly cou¬ 
pled A/” = 4 SYM plasma for T ~ 2 — 3Tc, which may re¬ 
sults in somewhat large Pj/s-y as shown in (13). However, 
OclT could vary drastically based on distinct models and 
approaches [50]. For example, by utilizing « 1/(2ctc), 
we find t-,, « 5 — 50 fm at T = 200 MeV given that 
UcjT « O.Ol — O.l ranging from lattice simulations to 
perturbative calculations. Although such approximated 
values of the relaxation time are smaller compared with 
the time scale of quark epoch in cosmology, it is inconclu¬ 
sive to pin down the exact values of rjj/s-y and in the 
cosmic plasma. On the other hand, our approach is not 
subject to the thermalized photons/leptons in the decon- 
fined phase. By knowing the photon-emission rate in the 
hadron gas, one could further analyze the photon/lepton 
transport in the nuclear matter. In fact, as shown in 
the phenomenological study in [JD], where the DG con- 
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ductivity of the interacting hadron gas is enhanced at 
low temperature, which may imply further suppression 
of lepton shear viscosity in the hadron epoch in cosmol¬ 
ogy. Nonetheless, the change of photon shear viscosity 
further requires the full photon spectrum in the hadron 
gas. The direct or indirect influence of the photon/lepton 
shear viscosity in the cosmic plasma on experimental ob¬ 
servables should be pursued in the future. 

Furthermore, on the theoretic side, one may incorpo¬ 
rate different effects such as a strong magnetic field and 
pressure anisotropy of the QGP, which may modify the 
conductivity of the QGP and the production rate of ther¬ 
mal photons laiiiHis] in various conditions. In addi¬ 
tion, for simplicity, we neglected chemical potentials and 
bulk fluctuations in our study. One could further general¬ 
ize the approach to include nonzero chemical potentials 
or investigate the bulk viscosity of photons/leptons in 
non-conformal plasmas. For example, as shown in the 
holographic model of QGD [IS], in which the photon- 
production rate is amplified by baryon chemical poten¬ 
tials particularly near the critical temperature, we may 
expect the shear viscosity of photons/leptons could be 
reduced near in the presence of chemical potentials. 
Also, although the imaginary part of the EM-current cor¬ 
relator in the sQGP is sufficient for evaluating the pro¬ 
duction rate and 77 .^, the real part also contributes to 


the self-energy of photons. The modification of the pho¬ 
ton dispersion relation in the strongly coupled medium 
may result in a negative refractive index The 

connection between the refractive index and shear vis¬ 
cosity of photons at strong coupling could make further 
impacts on different strongly interacting systems. More¬ 
over, the approach for analyzing the interplay between a 
perturabtive sector and a non-perturbative sector in this 
paper can also be applied to the semi-holographic model 
of QGD recently proposed in [50] . 
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